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Abstract 



This is an overview article on (/-special functions, a slightly extended version of an article 
to appear in the Encyclopedia of Mathematical Physics, Elsevier, 2006. 

1 Introduction 

In this article I give a brief introduction to g-special functions, i.e. (/-analogues of the classical 
special functions. Here q is a deformation parameter, usually < q < 1, where q = 1 is the 
classical case. The deformation is such that the calculus simultaneously deforms to a (/-calculus 
involving (/-derivatives and (/-integrals. The main topics to be treated are (/-hypergeometric 
series, with some selected evaluation and transformation formulas, and some (/-hypergeometric 
orthogonal polynomials, most notably the Askey- Wilson polynomials. In several variables we 
discuss Macdonald polynomials associated with root systems, with most emphasis on the A n 
case. The rather new theory of elliptic hypergeometric series gets some attention. While much 
of the theory of (/-special functions keeps q fixed, some of the deeper aspects with number 
theoretic and combinatorial flavor emphasize expansion in q. Finally we indicate applications 
and interpretations in quantum groups, Chevalley groups, afflne Lie algebras, combinatorics and 
statistical mechanics. 

As general references can be recommended [7j and |14j . 

Conventions 

q E C\{1} in general, but < q < 1 in all infinite sums and products. 
n, m, N will be nonnegative integers unless mentioned otherwise. 

2 g-Hypergeometric series 

Standard reference here is Gasper & Rahman 
2.1 Definitions 

For a, q € C the q-shifted factorial (a; q)k is defined as a product of k factors: 

(a-q) k :=(l-a)(l-aq)...(l-aq k - 1 ) {k G Z >0 ); (o;g) :=l. (2.1) 
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If \q\ < 1 this definition remains meaningful for k = oo as a convergent infinite product: 



(a; 9)00 := ]J(1 
j=0 



aq J 



We also write (01, . . . , a r ; for the product of r (/-shifted factorials: 

(d, . . . , a r ; q) k := (ai;q) k . • • (a r ; <?)fc (& G Z> or A; = 00). 



(2.2) 



(2.3) 



A q-hypergeometric series is a power series (for the moment still formal) in one complex vari- 
able z with power series coefficients which depend, apart from q, on r complex upper parameters 
01, . . . ,a r and s complex lower parameters b\, . . . , b s as follows: 



ai, 



,(f>s(ai, ■ ■ ■ ,a r ;bi, . . . ,b s ;q,z) 



fc=0 



(ai, . . . , a r ;q)k 
(61, . . .,b s ;q) k (q;q) k 



s—r+l 



(r,8€ Z> ). (2.4) 



Clearly the above expression is symmetric in ai, . . . , a r and symmetric in b\, . . . ,b s . On the 
right-hand side of (|2.4|) we have that 



(A; + l)th term (1 — aig* 



fl — a r 



,fc\s— r+l 



/cth term 



(1 - 61^) . . . (1 - M fe ) (1 - q k + l ) 



(2.5) 



is rational in q k . Conversely, any rational function in q k can be written in the form of the 
right-hand side of (|2.5|) . Hence, any series ^2^=0 c k W1 th cq = 1 and Ck+i/ck rational in g fc is of 
the form of a g-hypergeometric series (j2.4p . 

In order to avoid singularities in the terms of (|2.4|) we assume that 61, . . . , b s 1, q^ 1 ,q~ 2 , .... 
If, for some i, ai = q~ n then all terms in the series (|2.4|) with k > n will vanish. If none of the 
cij is equal to g _n and if |g| < 1 then the radius of convergence of the power series ()2.4j) equals 
00 if r < s + 1, 1 if r = s + 1, and if r > s + 1. 

We can view the g-shifted factorial as a g-analogue of the shifted factorial (or Pochhammer 
symbol) by the limit formula 



lim 



(g a ;q)k 

q"l (l-q) k 
Hence the q-binomial coefficient 



(a)k '■= a(a + 1) . . . (a + k — 1) 



{q;q)r 



{q;q)k(q;q)n-k 



(n,keZ, n > k > 0) 



tends to the binomial coefficient for q — * 1: 



lim 



77 






A 


q 


(I) 



(2.6) 



(2.7) 



(2.f 
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and a suitably renormalized (/-hypergeometric series tends (at least formally) to a hypergeometric 
series as q f 1: 



lim 



r+r' Ys+s' 



q ai ,...,q ar ,c 1 ,...,c r > 
Q bl ,- ■ ■,Q bs ,d 1 ,. ..,d s > 



; q, (q - 1 



^^.....a^ci-l)...^-!)^ 



fei, . . . , 6 S ' - 1) ... (4' - 1 
At least formally, there are limit relations between (/-hypergeometric series with neighbouring 



lim r (j) s 

a r — >oo 

lim r (j) s 
b s — >oo 







bi,. 


■ ,b s 


Oi, . . 






■ ,b s 



a r 

;q,b s z 



r -lC 



r<Ps-l 



ai, . . . , a r _i 
h,...,b s 
ai, . . . , a r 



y s -i 



(2.10) 
(2.11) 



A terminating (/-hypergeometric series Ylk=o c & zk rewritten as z n Ylk=o c n-kZ k yields an- 
other terminating (/-hypergeometric series, for instance: 



q ",Oi, ...,a s 



;q,z 



(— l) n q-\ n ( n+1 ) )—l ' ' ' ' a "' z n 

,b s ;q) n 

q~ n , q- n +\\ q^+% 1 q"+% ...b s 



q~ n+l a 1 1 , . . . ,q~ n+1 a s 1 ' ai . . . a s z 



(2.12) 



Often, in physics and quantum groups related literature, the following notation is used for 
(/-number, (/-factorial and (/-Pochhammer symbol: 



_ qz a - q 2 a 
qi — q 2 



fe fc-l 

II[7] ff> (W ff )fc:=n[° + ^ 

i=i j=0 



e z> ). 



(2.13) 



For (/ — ► 1 these symbols tend to their classical counterparts without the need for renormaliza- 
tion. They are expressed in terms of the standard notation (|2.1|) as follows: 



[*],! = ff-**<*-D -Sj^, ([a] q ) k = q'hHa-l) g -i*(*-i) 

2.2 Special cases 

For s = r — 1 formula (|2.4|) simplifies to 



(2.14) 



rVV-l 



ai, . . . , a r 
bi, . . . , 6 r _i 



fc=0 



(ai, . . . , a r ; g)fc 

(fei, . . .,b r -r,q)k (q;q)k 



(2.15) 



which has radius of convergence 1 in the non-terminating case. The case r = 2 of (|2.15|) is the 
(/-analogue of the Gauss hypergeometric series. 
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(/-Binomial series 



i9o{a;-;q,z) 



£ 

fc=0 



(a;q) k z k (az;q) c 



(q;q)k (z;q) c 



(if series is not terminating then \z\ < 1). (2.16) 



g-Exponential series 

e q (z) :=i^o(0; -;q,z) =J2j- 



OO h 



k=0 



{q;q)k (z;g) c 



\z\<D 



E q (z) := o 0o(-; -; q, -z) = ^ 



k=0 



(q;q)h 



e q {z) := 1 (j) 1 (0;-q2-q2 j -z) = ^ 



00 ifc(fc-l) 



/c=0 



(q;qh 



(-z;q)ao = (eq(-z)) 1 (z E C), 



[z G C). 



(2.17) 
(2.18) 
(2.19) 



Jackson's g-Bessel functions 

4 1] (x;q) : J qV+ \ q) °° (i x y 2 fc 
wj q)oo 

(<?;<?)oo 



0,0 _! 2 



■0 _I "+1_2 



(0 < x < 2), 



(2.20) 



u 1 2 

7 „+i ; ^ 



•2^;g)oc 4 1} (^9) (*>0), (2.21) 
(x > 0). (2.22) 



See (|3.26|) for the orthogonality relation for J„(x;q). 

If exp ? (z) denotes one of the three g-exponentials (|2.17j) - ()2.19j) then i(exp (J (ix) + exp (? (— ix)) 
is a ^-analogue of the cosine and — ^i(exp g (ix) — exp g (— ix)) is a ^-analogue of the sine. The 
three g-cosines are essentially the case v = — ^ of the corresponding g-Bessel functions (|2.2U|) — 
(I2.22|) . and the three g-sines are essentially the case v = \ of x times the corresponding g-Bessel 
functions. See also Suslov 



2.3 g-Derivative and g-integral 

The q- derivative of a function / given on a subset of R or C is defined by 

f{x) - f(qx) 



(DJ)(x) := 



(1 — g)x 



where x and gx should be in the domain of /. By continuity we set (D q f)(0) 
f'(0) exisits. If / is differentiable on an open interval / then 

lun(D q f)(x) = f'(x) (xel). 
?|1 



(2.23) 
/'(0), provided 

(2.24) 
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For a € R\{0} and a function / given on (0, a] or [a,0), we define the q-integral by 

oo oo 

f(x) d q x := a(l -q)J2 /(<"Z*) 9* = £ ^ ~ ( 2 " 25 ) 



fc=0 fc=0 



provided the infinite sum converges absolutely (for instance if / is bounded). If F(a) is given by 
the left-hand side of (|2.25|) then D q F = f. The right-hand side of Q2.25JI is an infinite Riemann 
sum. For q f 1 it converges, at least formally, to f{x) dx. 
For nonzero a, b G R we define 

b f(x) d q x := f f(x) d q x - f 1 f(x) d q x. (2.26) 
Jo Jo 

For a (/-integral over (0, oo) we have to specify a (/-lattice {aq k }u& for some a > (up to 
multiplication by an integer power of q) : 

/•a.oo 00 /•q~ n a 

/ f{x)d q x:=a{l-q) V /(ag fc ) g fc = lim / /(x) d,x. (2.27) 
Jo ^ "~Wo 

2.4 The g-gamma and g-beta functions 

The q- gamma function is defined by 

r q (z): J q '^^; q)1 ~ Z (^0,-1,-2,...) (2.28) 
[Q ,Q)oo 

t z ~ l E q (-(l- q)qt)d q t (Rez>0). (2.29) 



Then 



'V - ') = \z^ v M-. < 2 -' !( » 



limr„(z) = r(z). (2.32) 
The q-beta function is defined by 

t b-i (g*i g)oo , (Re6>0, o^0,-l,-2,...). (2-34) 
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2.5 The g-Gauss hypergeometric series 
q- Analogue of Euler's integral representation 

r 5 (c) f 1 +b-i ( t( i;q)oo (tzq a -,q) c 



9i{q a ,q b \q c ;q,z) 



T q (a)T q (c-b) J {tq c - b ;q)oo tz;q), 
By substitution of (|2.25|) . formula (|2.35f) becomes a transformation formula: 
20i (a, 6; c; g, z) = ; — ; 2 0i (c/o, z; az; g, 6). 

v^j q)oo { c 'i q)oo 

Note the mixing of argument z and parameters a, b, c on the right-hand side. 



d g t (Re b > 0, |z| < 1). 



Evaluation formulas in special points 

2 <f>i(a,b;c;q,c/(ab)) = 
2 Mq~ n ,b;c;q,cq n /b) = 
2 (f>i(q~ n ,b;c; q,q) = 

Two general transformation formulas 



(c/a,c/b; q) 
(c, c/(ab);q) c 
{c/b; q) n 



{\c/{ab)\ < 1), 



(c;q) n 

(c/W 
(c; g)„ 



291 



a, b 



{az;q)c 

{z;q)c* 

(abz/c; q) 
(z;q)oo 



292 



a, c/b 



c, az 



;q,bz 
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c/a,c/b abz 

;q, — 

c c 



Transformation formulas in the terminating case 
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q- n ,b 



392 



{c/b;q) r 
(c; q)n 

(q~ n bc~ 1 z;q) n 3 (p 



q n ,b,q n bc l z 
qt-Hc- 1 ^ ' ,<1,q 
q~ n , cb' 1 ^ 



c, qcb l z 1 ' ' 



(c/b; q) n 

s 301 

(c; q)n 



q n ,b,qz 

ql-nfc-l 



-1 



Second order (/-difference equation 

1 -q c 



z{q c -q a+b+l z)(D 2 u)(z) + 



l-q 



q"^r^ + q a K^—)z)(D q u)(z) 



l-q 



1 



1 - q a 1 



(2.35) 
(2.36) 



(2.37) 
(2.38) 
(2.39) 

(2.40) 
(2.41) 

(2.42) 
(2.43) 
(2.44) 



l-q 



u(z) = 0. (2.45) 



6 



Some special solutions of (|2.45j) are: 

2<M<? 



UlO) 

u 2 (z) 



a ,q b ;q c ;q,z), 



~1— c JL f.l+a— c „l+6— c. Jl—c. 

z 2<pi{q ,q ,q ,> 



z- a 2 Mq a ,q' 



a a— c+1 . „a— 6+1. 



9,9 



-a-6+c+l — 1 



(2.46) 
(2.47) 
(2.48) 



They are related by: 



(g^g^g^g)oc ( g b - 1 z,q 2 - b z- 1 ] g) 00 



(g fe - c z,g c - fe + 1 z- 1 ;g), 

(? l_c g a-6+l. g)oo ?c -a-6+l z -l. ^ z a 



{q l - b ,q a - c+1 ;q)oo (q b ~ c z,q 



b—c.y n c— 6+1 ~— 1. 



u 3 (z). (2.49) 



2.6 Summation and transformation formulas for 



! series 



An r 4>r-i series (j2.15|) is called balanced if b\ . . . 6 r _i = gai . . . a r and z = q, and the series is 

i 

called very-well-poised if oai = 0261 = 0362 = ■ ■ ■ = a r 6 r -i and qal = 02 = —03. The following 
more compact notation is used for very-well-poised series: 

r W r _i(ai;a 4 ,a 5 , . . . ,a r ;q,z) := r <t>r-i 



1 1 
ai,qa£,-qaZ,a4, 



II ;9,£ 
af , -a{ , 901/04, . . . , gai/a r 



(2.50) 



Below only a few of the most important identities are given. See [7J for many more. An important 
tool for obtaining complicated identities from more simple ones is Bailey's Lemma, which can 
moreover be iterated (Bailey chain), see Ch.3]. 



The g-Saalschiitz sum for a terminating balanced 



3<P2 



392 



a,b,q n 



c, q 1 n abc 1 ' 



(c/a,c/b; q) n 
(c, c/(ab);q) r 



(2.51) 



Jackson's sum for a terminating balanced %W>j 



W 7 (a;b,c,d,q n+1 a 2 /(bcd),q- n ;q,q) 



(qa, qa/ (be), qa/ (bd),qa/ (cd); q) r 



(qa/b, qa/c, qa/d, qa/ (bcd);q) n 
Watson's transformation of a terminating %W-j into a terminating balanced 4 </>3 



(2.52) 



^n+2„2 



iW 7 a;b,c,d,e,q n ;q, 



a 2 \ (qa,qa/(de);q) 7 



bede 



(qa/d, qa/e;q) n 



4<?>3 



q n ,d,e,qa/(bc) 
qa/b, qa/c, q~ n de/a 



5 9,9 



. (2.53) 



7 



Sears' transformation of a terminating balanced 4^3 



4<?3 



q n ,a,b,c^ 
d,e,f 



(e/a,f/a;q) n 

— / — 7-^ a 4<P3 

{e,f;q)n 



q n ,a,d/b,d/c 
d, q^a/e^^a/f' 



(2.54) 



By iteration and by symmetries in the upper and in the lower parameters, many other versions 
of this identity can be found. An elegant comprehensive formulation of all these versions is as 
follows. 

Let X1X2X3X4X5X6 = q l ~ n . Then the following expression is symmetric in x\, x 2 , X3, x±, x§, xq: 



X1X2X3X4, X 1 X 2 XsX 5 , X1X2X3XQ; 



q '%x 2 x 3 ,x l x 3 ,xix 2 . . 

403 ;q,q • (2.55) 

(XiX2X3) n \_XiX 2 X 3 Xi,XiX 2 X 3 X 5 ,XiX 2 X 3 X & 

Similar formulations involving symmetry groups can be given for other transformations, see |27| . 

Bailey's transformation of a terminating balanced 10W9 

(qa, qaj (ef), (qa) 2 / (bcde), (qa) 2 / (bcdf); q) n 



q n+2 a 3 

10 W 9 I a; b, c, d, e, /, , , , , q n ; q, 



bcdef ' 



(qa/e,qa/f, (qa) 2 /(bcdef), (qa) 2 /(bcd);q) n 



x io W 9 



n+2 a 3 



qa qa qa qa 
bed ' cd ' bd ' be ' ' ' bcdef 



<q ;q,q 



2.7 Roger s-Ramanuj an identities 



OO u2 

o0i(-;O;?,g) = 



(q;q)h (<?,<? 4 ;<? 5 )c 



o<M-;0;<?,<? 2 ) = ^2 



k=0 

00 q k(k+l) 



k=0 



(q;q)k (<?W;<? 5 )c 



2.8 Bilateral series 
Definition (|2.1|) can be extended by 



(a;q) k ■-- 



(a;q)c 



(k G Z) 



(aq k ;q)oo 

Define a bilateral q-hypergeometric series by the Laurent series 



-IP* 



ai,...,a r 
bi,...,b s 

CO 



\q,z 



rip s (ai, ...,a r ;bi,.. .,b s ;q,z) 



k=—oo 



(61, . . .,b s ;q) k 



(2.56) 



(2.57) 
(2.58) 



(2.59) 



E ( ° 1 '---'^ ;g , )fc ((-l)V fc(fe - 1) )° V (a 1 ,...,a r ,b 1 ,...,b s ^0,s>r). (2.60) 



The Laurent series is convergent if \b± . . . b s j (a\ . . . a r )\ < \z\ and moreover, for s = r, \z\ < 1. 
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Ramanajan's yipi summation formula 



(q,c/b,bz,q (bz^q)^ 

itpi(b c;q,z) = - — ( c/6 < z < 1). (2.61) 

(c, g/6,z,c/(6z);g) 00 

This has as a limit case 

0^(-;c; g ,*)= (g /^ ; . g)o ° (N>|c|), (2.62) 

(C, C/ZJ^oo 

and as a further specialization the Jacobi triple product identity 

oo 

ST (-l) k ga*(fc-i) / = ( g> g)oo ( z ^ o), (2.63) 

fc=— oo 

which can be rewritten as a product formula for a i/ieta function: 

oo oo 

# 4 ( x; g ) := ("l) fc ? fe2 e 2 ™ fcx = II (1 - 9 2fc ) (1 - 2g fc ~ 1 cos(27rx) + q 4k ~ 2 ) . (2.64) 

fc=— oo fc=l 

3 g-Hypergeometric orthogonal polynomials 

Here we discuss families of orthogonal polynomials {p n (x)} which are expressible as terminating 
g-hypergeometric series (0 < q < 1) and for which either (i) P n {x) := p n {x) or (ii) P n {x) := 
Pn(^{x + x~ 1 )) are eigenfunctions of a second order ^-difference operator, i.e.: 

A(x) P n {qx) + B(x) P n (x) + C(x) P n (q- l x) = X n P n (x), (3.1) 

where A(x), B(x) and C{x) are independent of n, and where the A n are the eigenvalues. The 
generic cases are the four-parameter classes of Askey-Wilson polynomials (continuous weight 
function) and q-Racah polynomials (discrete weights on finitely many points). They are of type 
(ii) (quadratic g-lattice). All other cases can be obtained from the generic cases by specialization 
or limit transition. In particular, one thus obtains the generic three-parameter classes of type 
(i) (linear ^-lattice). These are the big q- Jacobi polynomials (orthogonality by g-integral) and 
the q-Hahn polynomials (discrete weights on finitely many points). For all these families the 
standard formulas are given in Koekoek & Swarttouw [TP] . 

3.1 Askey-Wilson polynomials 

These were introduced by Askey &: Wilson [I] , 



Definition as g-hypergeometric series 

( a\ r a h j\ \ {ab,ac,ad;q) n 

p n (cos6) = Pn{cos 0; a, b, c,d q) := 4 </> 3 

a n 

This is symmetric in a, b, c, d. 



q n ,q n 1 abcd,ae ld ,ae %e 
ab, ac, ad 



(3.2) 
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Orthogonality relation Assume that a, b, c, d are four reals, or two reals and one pair of 
complex conjugates, or two pairs of complex conjugates. Also assume that \ab\, \ac\, \ad\, \bc\, \bd\, 
\cd\ < 1. Then 

/ Pn{x) Pm(x) w(x) dx + y^Pn(^fc) Pmjxk) ^k = K 8 n>m , (3.3) 

k 

(e™;q) 



where 



2tt sinOw (cos 9) 



IOC: 



{ae id ,be id ,ce ie ,de ie ;q) c 
(abed; q)^ h n 1 — abedq 71 ^ 1 (q, ab, ac, ad, be, bd, cd; q) r 



(3.4) 
(3.5) 

u (q, ab,ac, ad, bc,bd,cd; q)^ ' ho 1 — abcdq 2n ~ l (abcd;q) n 

and the x k are the points \(eq k + e~ 1 q~ k ) with e any of the a, b, c, d of absolute value > 1; the 
sum is over the k £ Z>o with \eq k \ > 1. The Uf. are certain weights which can be given explicitly. 
The sum in (|3.3|) does not occur if moreover \a\, \b\, |c|, \d\ < 1. 

A more uniform way of writing the orthogonality relation 1)3.3(1 is by the contour integral 

— ^PnUz + Z )) Pm (- 2 (Z + Z )) {aZiaz - lMbz - 1 ^- 1} ^- 1 . q)oo ~ = ^nS n , m , 

(3.6) 

where C is the unit circle traversed in positive direction with suitable deformations to separate 
the sequences of poles converging to zero from the sequences of poles diverging to oo. 
The case n = m = of 1)3.6)) or 1)3.3)1 is known as the Askey-Wilson integral. 

(/-Difference equation 

A(z)P n (qz) - (A(z) + Aiz-^Pniz) + A(z" 1 )P n (q~ 1 z) = (q~ n - 1)(1 - q n - 1 abcd)P n (z), (3.7) 
where P n {z) = Pn (\(z + z" 1 )) and A(z) = (1 - az)(l - bz)(l - cz)(l - dz)/((l - z 2 )(l - qz 2 )) . 

Special cases These include the continuous q-Jacobi polynomials (2 parameters), the contin- 
uous q-ultraspherical polynomials (symmetric one-parameter case of continuous g-Jacobi), the 
Al-Salam-Chihara polynomials (Askey-Wilson with c = d = 0), and the continuous q-Hermite 
polynomials (Askey-Wilson with a = b = c = d = 0). 

3.2 Continuous g-ultraspherical polynomials 

Definitions as finite Fourier series and as special Askey-Wilson polynomial 

C n (cos#;/3 | q) := y (#g)fc(#g)"-fc e i(n-2k)0 (3 8) 

W,q) k (q;q) n - k 

= ^^Pn(cos^MM,-M,-^M I ?)• (3-9) 

Wi Q)n 
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Orthogonality relation (—1 < (3 < 1) 



— / ( '„(cos#: i.q)( ',„(cos#: i.q) 



(f^qthqU 1-/3 (/3 2 ;g)n x 
aa — -r^. ; — tt— — ; — o n ,m- 



(P 2 ,q;q)oo i - Pq n {q;q) n 



(3.10) 



g-Difference equation 

A(z)P n (qz) - (A(z) + A(z~ 1 ))P n (z) + A{ Z - l )P n {q- 1 z) = (q~ n - 1)(1 - q n f3 2 )P n {z), (3.11) 
where P n (z) = C n {\{z + z' 1 ); (5 \ q) and A(z) = (1 - /?z 2 )(l - q(3z 2 )/((l - z 2 )(l - qz 2 )). 

Generating function 

(Pe t9 z,pe- ie j 1 q)^ = j, c ^ cose _ f3 ^^ zn (| z | < ^ < < tt, -1< /3 < 1). (3.12) 

n=0 



(e if? z,e- ie z; 



Special case: the continuous g-Hermite polynomials 

H n (x | q) = (q;q) n C n (x;0 \ q). 



Special cases: the Chebyshev polynomials 



C n (cos9;q \ q) = U n (cos6) :- 



sin((n+ 1)0) 
sin 6 



lim ^ Q \ n C„(cos 0;P\q) = T„(cos 6) := cos(n#) [n > 0). 
3.3 g-Racah polynomials 

Definition as g-hypergeometric series (n = 0, 1, . . . , N) 

' n ,al3q n+l ,q-y,-f5qy +l 



(3.13) 

(3.14) 
(3.15) 



Rn(q- y +lSq y+1 ; a, £7, 5 | g) := 4</>3 



Orthogonality relation 



;q,q 



(a, (35 ot j = q 



_ r ,-N-i\ 



(3.16) 



^ i?„( 9 -^ + ^q V+1 ) R m {q~ V + 7<V +1 ) W„ = Mn.m, 
j/=0 



(3.17) 



where a;^ and /i n can be explicitly given. 
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3.4 Big g-Jacobi polynomials 
Definition as g-hypergeometric series 

P n (x) = P n {x; a, b, c; q) := 3 (/> 2 



q n ,q n+1 ab,x 
qa, qc 



(3.18) 



Orthogonality relation 

r P n (x) P m (x) ( ° l ,C _ x X] g) °° d q x = h n S n>m , (0<a<q-\ 0<b<q-\ c<0), (3.19) 
where h n can be explicitly given. 

(/-Difference equation 

A(x)P n (qx) - (A(x) + C(x))P n (x) + C(x)P n (q- 1 x) = (q~ n - 1)(1 - abq n+1 ) P n {x) , (3.20) 
where A{x) = aq(x — l)(bx — c)/x 2 and C{x) = (x — qa)(x — qc)/x 2 . 

Limit case: Jacobi polynomials P^^\x) 

lim P n (x; q a , q p , -q~ l d; q) 



<?Ti {a + 1)„ 

Special case: the little g-Jacobi polynomials 

P 



nl p(a,/3) ' 2X + d ~ 1 
n 



d+ 1 



(3.21) 



(x; a, b; q) = (-6)-"g-^ (n+1) p n ( gte; ft, a , 0; , 

{qa; q)n 

= 2<Pi{q~ n ,q n+1 ab;qa;q,qx), 



(3.22) 
(3.23) 



which satisfy orthogonality relation (for < a < q and b < q ) 



p n {x; a,b; q) p m (x; a,b; q) , 
n g) 



(gaj; g)oo w a , 
x B 9 a„x 



{q,qab;q)oo {l-q){qa) n (q,qb;q) n 



(ga, g6; g)^ 1 - abq 2n+1 {qa, qab; g) r 



(3.24) 



Limit case: Jackson's third g-Bessel function (see ()2.22(l and jTH]) 



hm PN-n{q N+k ;q u ,b;q) = %? )o ° g^"+ fc ) j( 3 >(2g^+ fc ); g) {v > -1), (3.25) 



iV— ►oo 



by which Q3.24J1 tends to the orthogonality relation for J„(x;q): 

oo 

J2 4 3) (2^ (n+fe) ;<7)4 3) (2g5 (m+fc) ;g)g fc = <5 n , m g- n (n,meZ). 



(3.26) 



k=— oo 



12 



3.5 g-Hahn polynomials 
Definition as g-hypergeometric series 

~q~ n ,q 



Q n (x;a,(3,N;q) := 3 2 



n » n+1 al3,x 



qa,q 



-N 



(n = 0,l,...,N). 



Orthogonality relation 



N 



y=0 

where h n can be explicitly given. 

3.6 Stieltjes-Wigert polynomials 
Definition as q-hypergeometric series 

S n (x;q) 



(q;q)r 



o 



;q,-q n+1 x 



The orthogonality measure is not uniquely determined: 

1 



S n (q2 X ;q) S m (q*x;q)w(x)dx 



l (q;q)r 



J n,mi 



where, for instance, 



w(x) 



i 

q-2 



or 



q 2 



exp 



log 2 x 
"21og(g-i; 



log^" 1 ) (<?, -qix^-qix-^q)^ y/2ir\og{q x ) 

3.7 Rahman- Wilson biorthogonal rational functions 

The following functions are rational in their first argument: 

R n (h(z + z^ 1 ); a, b, c, d,e) := ioWg(a/e; qj (be),q/ (ce),q/ (de), az, a/z, q n ~ 1 abcd 
They satisfy the biorthogonality relation 

Rn (\{z + z~ l );a,b,c, d,e) R m (\(z + z~ l ); a, b, c, d, 

where the contour C is as in (j3,6|) . and where 

(z 2 , z~ 2 , abcdez, abode/ z; q) 



j . dz 
w(z ) — = 2h 
abode J z 



w(z) 



(az, a/z, bz, b/z, cz, c/z, dz, d/z, ez, e/z; q)oo ' 

(bcde, acde, abde, abce, abed; q)oo 
(q, ab, ac, ad, ae, be, bd, be, cd, ce, de; q)^ ' 
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and h n /h(j can also be given explicitly. For ab = q N , n,m £ {0, 1, . . . , N} there is a related 
discrete biorthogonality of the form 

N 

Rn (\{aq k + a' 1 q~ k );a, b, c, d, ej R m (±{aq k + a' l q~ k )\ a, b, c, d, w k = (n^m). 

(3.35) 

See Rahman [23| and Wilson |2S|. 



4 Identities and functions associated with root systems 
4.1 ^-Function identities 

Let R be a root system on a Euclidean space of dimension /. Then Macdonald ^H] generalizes 
Weyl's denominator formula to the case of an affine root system. The resulting formula can be 
written as an explicit expansion in powers of q of 

oo 
n=l 

which expansion takes the form of a sum over a lattice related to the root system. For root system 
Ai this reduces to Jacobi's triple product identity (|2.63|) . Macdonald's formula implies a similar 
expansion in powers of q of rj(q) l+ \ R \, where 77(g) is Dedekind's rj-function rj(q) := <? 1//24 (g; q)oo- 



4.2 Constant term identitities 



Let R be a reduced root system, R + the positive roots and k € Z>o- Macdonald ^Hl conjectured 
the second equality in 



IrUa^R+i 6 a ;q)k(qe a ;q)kdx 

jrp dX 



I r 



cr{ n rid 

a£R+ i=l 



e- a )(l-q l e a ) 



n 

i=i 



kdi 
k 



(4.1) 



where T is a torus determined by R, CT means the constant term in the Laurent expansion 
in e a , and the d{ are the degrees of the fundamental invariants of the Weyl group of R. The 
conjecture was extended for real k > 0, for several parameters k (one for each root length), and 
for root system BC n , where Gustafson's |Hj 5-parameter n- variable analogue of the Askey- Wilson 
integral f (|3,6j) for n = 0) settles: 



[0,27r] n 



|A(e* e \...,e^)| 2 ^^ 



n 

2 n n! J] 

i= 



(i, t n+ i- 2 abcd; q) c 



where 



A(*):= II 



ZiZj, Zij Zj] q)oo j-j- ( z |;?)oo 
\tZiZj 5 tZij Zjj £/)oo ( CI Z j , C2^j, dZj } q)oo 



(4.2) 



(4.3) 



l<i<j'<n x J ' "' J ' *' j = 

Further extensions were in Macdonald's conjectures for the quadratic norms of Macdonald poly- 
nomials associated with root systems (see ^4.4|) . and finally proved by Cherednik [§]. 
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4.3 Macdonald polynomials for root system Ai-i 

Reference for this section is Macdonald |2U1 Ch. VI]. Let n € Z>o- We work with partitions 
A = (Ai, . . . , X n ) of length < n, where Ai > • • • > A n > are integer. On the set of such partitions 

we take the partial order A < fi -4=^- Ai + • • • + A n = \i\ H V fi n and Ai + • • • + Aj < fj,\ + • • • + //j 

(i = 1, ... ,n — 1). Write A < [i iff A < fi and A / /i. The monomials are z a = z" 1 . . . z" n 
(ai, . . . , a n € Z>o). For A a partition the symmetrized monomials m\(z) and the Schur functions 
s\(z) are defined by: 



) := zOL (sum over all distinct permutations a of (Ai, . . . , A n )), (4-4) 



a 

(z) det(^ + " J )^ J= i,...,n (4 } 

det(^ )ij=i,..., n 

We integrate a function over the torus T := {2 G C™ | = . . . = \z n \ = 1} as 

/(*) := ^ ... f(e ie \. . . , e^) d6 1 . . . dd n . (4.6) 

Definition For A a partition and for < t < 1 the (analytically defined) Macdonald polynomial 
P\(z) = P\(z; q, t) is of the form 

P\(z) = P x (z; q, t) = m\(z) + u\^m^(z) (u AjjU G C) 



such that for all // < A 



J P x (z)m fl (z)A(z)dz = 0, 



where 



Orthogonality relation 

(q x i- x n^ i ,q Xi -^ +1 t j - i ;q) c 



1 z 1 / n Ai— Ajj-J— % n Ai—Aj-\-l-i-j—l. n \ 



g-Difference equation 



EII^v^;^) = (X> Al ^ J (4-9) 



i=l J \i=± 



where T qtZi is the q-shift operator. T q>Zi f(zi, . . . , z n ) := f(zi, . . . , qzi, . . . , z n ). See (2U1 Ch. VI, 
§3] for the full system of g-difference equations. 
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Special value 

p a(1 , *, . ..,«->, q , t) = n n ■ < 4 - 10 ) 

i=l i<j W ' ( /A ! -A J 

Restriction of number of variables 

f , Ai,A 2 ,...,A n _i,ci0i ! • • • ,z n -i,0;?,i) = Pai,a 2 ,...,A„-i(2i ; • • ■ ,2 n _i;g,i). (4.11) 

Homogeneity 

i 3 A 1 ,...,A„(^;g,0 =^l---- z n^A 1 -l,...,A rl -l(-2;g,t) (A n >0). (4.12) 



Self-duality Let A, [i be partitions. 

P x (q^H n - 1 ,q^t n - 2 , ...,q^;q,t) _ P^(q Xl t n -\ q x H n - 2 , . . . ,q Xn ;q,t) 



(4.13) 



Special cases and limit relations 

Continuous q-ultraspherical polynomials (see (|3.8|0 : 

P m , n (re ld , re- 19 ; q, t) = ^ q)m ~ n r m+n C m _ n (cos 9;t\q). (4.14) 

V^i Q)m—n 

Symmetrized monomials (see (|4.4|) ): 

P x (z;q,l)=m x (z). (4.15) 

Schur functions (see (|4.5(l ): 

P x (z;q,q) = s x (z). (4.16) 
Hall-Littlewood polynomials (see |2U1 Ch. Ill]): 

Px(z;0,t) =P x {z;t). (4.17) 

Jack polynomials (see |2(J1 §VI.10]): 

limP x (z;q,q a )=P { x 1/a \z). (4.18) 

Algebraic definition of Macdonald polynomials Macdonald polynomials can also be de- 
fined algebraically. We work now with partitions A (Ai > A2 > • • • > 0) of arbitrary length 
Z(A), and with symmetric polynomials in arbitrarily many variables X\,X2, ■ ■ ■ , which can be 
canonically extended to symmetric functions in infinitely many variables xi,x%, . . . . The rth 
power sum p r and the symmetric functions p x are formally defined by 

Pr = ^ x h P\ =P\xP\2 ••• • ( 4 - 19 ) 
i>l 
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Put 



z\ := || i m% mil , where mi = mj(A) is the number of parts of A equal to i. (4.20) 

i>l 

Define an inner product ( , ) Q) t on the space of symmetric functions such that 

IW , x 

(P\,Pfl)q,t = $\,n z\ [[ _ fX . ■ (4.21) 

1=1 

For partitions \,fi the partial ordering A > \i means now that Ylj>i^j = Ylj>i^j an d + 

h A$ > jtii H \- fj,i for all i. The Macdonald polynomial P\(x; q,t) can now be algebraically 

defined as the unique symmetric function P\ of the form P\ = J2fi<\ u \,ti m n ( n A,/i £ C, it^.A = 1) 
such that 

(Pa,^)^ = if A y^. (4.22) 

If l(X) < n then the newly defined P\(x) with x n+ \ = x n+ 2 = . . . = coincides with P\(x;q,t) 
defined analytically, and the new inner product is a constant multiple (depending on n) of the 
old inner product. 



Bilinear sum 



Generalized Kostka numbers The Kostka numbers K\ t u occurring as expansion coefficients 
in s\ = Ylu^\fi m n were generalized by Macdonald to coefficients K\^(q,t) occurring in con- 
nection with Macdonald polynomials, see [2UJ §VI.8]. Macdonald's conjecture that K\ tlM (q,t) is 
a polynomial in q and t with coefficients in Z>q was fully proved in j^j. 



4.4 Macdonald-Koornwinder polynomials 

Macdonald |5J also introduced Macdonald polynomials associated with an arbitrary root sys- 
tem. For root system BC n this yields a three-parameter family which can be extended to 
the five-parameter Macdonald-Koornwinder polynomials (M-K polynomials) ^3]. They are or- 
thogonal with respect to the measure occurring in (|4.2|) with A(z) given by (|4.3|) . The M-K 
polynomials are re- variable analogues of the Askey- Wilson polynomials. All polynomials just 
discussed tend for q f 1 to Jacobi polynomials associated with root systems. 

Macdonald conjectured explicit expressions for the quadratic norms of the Macdonald poly- 
nomials associated with root systems and of the M-K polynomials. These were proved by 
Cherednik by considering these polynomials as Weyl group symmetrizations of non-invariant 
polynomials which are related to double affine Hecke algebras. See Macdonald 22 . 
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5 Elliptic hypergeometric series 

See Gasper & Rahman [7[ Ch. 11] and references given there. 

Let p, q € C, \p\, \q\ < 1. Define a modified Jacobi theta function by 

6(x;p) := (x,p/x;p)oo (i / 0), (5.1) 

and the elliptic shifted factorial by 

(a;q,p) k := 0(a;p)9(aq;p) ...9(aq k ~ 1 ;p) (keZ >0 ); (a;q,p) := 1, (5.2) 



(ai, . . .,a r ;q,p) k := {a 1 ;q,p) k . . . (a r ;q,p) k , (5.3) 
where a, ai, . . . , a r ^ 0. For q = e 27rW7 , p = e 27nT (Imr > 0) and a S C we have 



ae 



' ' 1, 1 ^ „ . = -a" 1 ^. (5.4) 



A series X]^lo Cfc w ith c k +\jc k being an elliptic (i.e. doubly periodic meromorphic) function 
of k considered as a complex variable, is called an elliptic hypergeometric series. In particular, 
define the r E r -i theta hypergeometric series as the formal series 

r E r _ x {a x ,.. ■ ,a r ;bi, . . . , b r -r, q,p; z) := > 77 — -. — . (5.5) 

^ (61, . . . ,b r -i;q,p) k {q;q,p) k 

It has g(k) := c k+ \/c k with 

z9{aiq x ]p) . . .6(a r q x ;p) 



x+1. 



p) 6(biq x ;p) . . . 0(b r ^iq x ;p) 



By (|5.4[) 5(3;) is an elliptic function with periods cr _1 and tct~ 1 (q = e 2ma , p = e 2mT ) if the 
balancing condition a% . . . a r = qb\ . . . b r -\ is satisfied. 

The r V r -\ very-well-poised theta hypergeometric series (a special r E r -\) is defined, in case 
of argument 1, as: 

T/ / \ ^0{aiq 2k ]p) (a 1 ,a 6 ,...,a r ;q,p) k q k 

rV r -. 1 [a 1 \ae,...,a r ;q,p) := > — —- r? -. (5.6) 

9{ai;p) {qa 1 /as,...,qai/a r ;q,p)k {q;q,P)k 

The series is called balanced if a| . . . a? r = a[~ 6 q r ~' i . The series terminates if, for instance, 
a r = q~ n . 

Elliptic analogue of Jackson's gWj summation formula ()2.52j) 

T7- / 7 , n+1 2/, h a -n \ (qa,qa/(bc),qa/{bd),qa/(cd);q,p) n 

w V 9 {a;b,c,d,q ^ a {bcd),q ;q,p) = — — — r — . (5.7 

[qa/b, qa/c, qa/d, qa {bed); q,p) n 
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Elliptic analogue of Bailey's ioWg transformation formula (|2.56j) 



12V11 (a-bcdef q ~ n - q p] = (g Q »W(e/), (qa) 2 / (bcde) , (ga) 2 /(bcdf); q,p) n 

12 11 V bcdef ' ' ' J (qa/e,qa/f,(qa) 2 /(bcdef),(qa) 2 /(bcd);q,p) n 

X 12V11 t—;; — -, ti, t-, e, /, , g g,p . (5.8) 

\oca ca oa oc bcdej J 

Suitable 12^11 functions satisfy a discrete biorthogonality relation which is an elliptic analogue 
of (HOST) . 

Ruijsenaars' elliptic gamma function (see Ruijsenaars [21] ) 

~ 1 _ z - 1 o^'+ 1 » fc + 1 

j,fc=0 

which is symmetric in p and q. Then 

F(qz;q,p) = 9(z;p)r(z;q,p), F(q n z;q,p) = (z;q,p) n T(z;q,p). (5.10) 



6 Applications 

6.1 Quantum groups 

A specific quantum group is usually a Hopf algebra which is a g-deformation of the Hopf algebra 
of functions on a specific Lie group or, dually, of a universal enveloping algebra (viewed as Hopf 
algebra) of a Lie algebra. The general philosophy is that representations of the Lie group or 
Lie algebra also deform to representations of the quantum group, and that special functions 
associated with the representations in the classical case deform to g-special functions associated 
with the representations in the quantum case. Sometimes this is straightforward, but often new 
subtle phenomena occur. 

The representation theoretic objects which may be explicitly written in terms of g-special 
functions include matrix elements of representations with respect to specific bases (in particular 
spherical elements), Clebsch-Gordan coefficients and Racah coefficients. Many one- variable q- 
hypergeometric functions have found interpretation in some way in connection with a quantum 
analogue of a three-dimensional Lie group (generically the Lie group SL(2, C) and its real forms). 
Classical by now are: little g-Jacobi polynomials interpreted as matrix elements of irreducible 
representations of SU q {2) with respect to the standard basis; Askey- Wilson polynomials similarly 
interpreted with respect to a certain basis not coming from a quantum subgroup; Jackson's third 
q-Bessel functions as matrix elements of irreducible representations of E q (2); q-Hahn polynomials 
and (/-Racah polynomials interpreted as Clebsch-Gordan coefficients and Racah coefficients, 
respectively, for SU q {2). See for instance Vilenkin & Klimyk [28 : and Koelink |llj 

Further developments include: Macdonald polynomials as spherical elements on quantum 
analogues of compact Riemannian symmetric spaces; q- analogues of Jacobi functions as matrix 
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elements of irreducble unitary representations of SU q (l, 1); Askey- Wilson polynomials as matrix 
elements of representations of the SU(2) dynamical quantum group; an interpretation of discrete 
12V11 biorthogonality relations on the elliptic U(2) quantum group (see [T2]^. 

Since the g-deformed Hopf algebras are usually presented by generators and relations, iden- 
tities for g-special functions involving non-commuting variables satisfying simple relations are 
important for further interpretations of g-special functions in quantum groups, for instance: 

g-Binomial formula with (/-commuting variables 

y n-k x k ^ xy _ qyx y -q 

Functional equations for q- exponentials with xy = qyx 

e q (x + y) = e q (y)e q (x), E q (x + y) = E q (x)E q (y), (6.2) 
e q (x + y-yx) = e q (x)e q (y), E q (x + y + yx) = E q (y)E q (x). (6.3) 

See |15| for further reading. 



6.2 Various algebraic settings 

Classical groups over finite fields (Chevalley groups) 

(/-Halm polynomials and various kinds of (/-Krawtchouk polynomials have interpretations as 
spherical and intertwining functions on classical groups (GL n , SO n , Sp n ) over a finite field F q 
with respect to suitable subgroups, see Stanton [2*5] . 



Affine Kac-Moody algebras (see [T7] ) 

The Rogers- Ramanuj an identities (|2.57[) . (|2.58j) and some of their generalizations were inter- 
preted in the context of characters of representations of the simplest affine Kac-Moody algebra 
. Macdonald's generalization of Weyl's denominator formula to affine root systems (see iJ4.3|) 
has an interpretation as an identity for the denominator of the character of a representation of 
an affine Kac-Moody algebra. 



6.3 Partitions of positve integers 

Let n be a positive integer, p(n) the number of partitions of n, pn(ti) the number of partitions 
of n into parts < N, pdist(^) the number of partitions of n into distinct parts, and p dd(^) the 
number of partitions of n into odd parts. Then Euler observed: 

1 00 00 

___ = J3p(n)g», ^^ = Y^p N {n)q\ (6.4) 

00 | :x. 

{-q\q)oo = ^Pd\st(n)q n , = ^2Podd{n)q n , (6.5) 

n=0 [q,q >QO n=0 
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and 

(-q;q)oo = — or— , Pdist(w) = p Q dd(w). (6.6) 
(?;r)oo 

The Rogers- Ramanuj an identity (|2,57|) has the following partition theoretic interpretation: 
The number of partitions of n with parts differing at least 2 equals the number of partitions of 
n into parts congruent to 1 or 4 (mod 5). Similarly, (j2,58|) yields: The number of partitions of 
n with parts larger than 1 and differing at least 2 equals the number of partitions of n into parts 
congruent to 2 or 3 (mod 5). 

The left-hand sides of the Rogers- Ramanuj an identities (|2.57[) and (|2.58|) have interpretations 
in the Hard Hexagon Model, see (Hj- Much further work has been done on Rogers- Ramanuj an type 
identities in connection with more general models in statistical mechanics. So-called fermionic 
expressions do occur. 
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